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Abstract 

Background: Alternative sigma factors trigger various adaptive responses. Lactobacillus sakei, a non-sporulating 
meat-borne bacterium, carries an alternative sigma factor seemingly orthologous to of Bacillus subtilis, best 
known for its contribution to the initiation of a large starvation response ultimately leading to sporulation. As the 
role of a^-like factors has been little studied in non-sporulating bacteria, we investigated the function of in L. 
sakel. 

Results: Transcription of sigti coding for was hardly affected by entry into stationary phase in our laboratory 
conditions. Twenty-five genes potentially regulated by in L sakei 23 K were revealed by genome-wide 
transcriptomic profiling of sigti overexpression and/or quantitative PCR analysis. IVlore than half of them are 
involved in the synthesis of a DNA uptake machinery linked to genetic competence, and in DNA metabolism; 
however, overproduction did not allow detectable genetic transformation, was found to be conserved in 
the L. sakei species. 

Conclusion: Our results are indicative of the existence of a genetic competence state activated by in L. sakel, 
and sustain the hypothesis that a'^-like factors in non sporulating Firmicutes share this common function with the 
well-known ComX of naturally transformable streptococci. 



Background 

Sigma factors are subunits of the RNA polymerase com- 
plex responsible for specific recognition and melting of 
promoter DNA, which enable the polymerase to initiate 
transcription. All eubacteria of known genome sequence 
code for at least one sigma factor, called primary, house- 
keeping or vegetative, and most encode additional sigma 
factors. For example, Streptomyces coelicolor or Somn- 
gium cellulosum carry as many as 60 to 80 predicted 
sigma factors [1,2]. These so-called alternative sigma 
factors may be induced or activated by specific environ- 
mental signals, and consequently redirect transcription 
by competitively associating with the core RNA poly- 
merase. Alternative sigma factors have been shown to 
mediate various cellular responses linked to stress con- 
ditions, growth transitions or morphological changes 
and development [1]. 

Sigma factors are classified into two structurally and 
evolutionarily distinct superfamilies [3], a 
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Most eubacterial sigma factors belong to the a super- 
family, which is further divided into four phylogenetic 
groups on the basis of protein structure and physiologi- 
cal function [1,4]. Group 1 comprises the housekeeping 
sigma factors. Group 2 is close to group 1 but accom- 
modates non essential sigma factors, including the mas- 
ter regulator of general stress response in stationary 
phase, RpoS, as was well characterized in Escherichia 
coli. Sigma factors in group 3 are phylogenetically 
diverse, and regulate major cellular functions such as 
sporulation, motility, heat-shock or general stress 
response. Group 4, known as the extracytoplasmic func- 
tion (ECF) subfamily, has been distinguished more 
recently. It comprises highly diverged sigma factors 
mainly involved in responses to extracytoplasmic sti- 
muli, which may affect the correct folding of envelope 
proteins. These factors typically contain only domains 
refgrped to as 2 and 4, involved in core polymerase 
binding and promoter DNA recognition and melting [3], 
with a spacer domain of less than 50 residues [2]. How- 
ever, due to the high divergence across sigma factors, 
their classification in the previously identified 
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phylogenetic groups may need to be revised, and new 
cellular functions controlled by sigma factors may be 
discovered [4]. 

Our research concerns a putative factor in the lac- 
tic acid bacterium Lactobacillus sakei. The closest char- 
acterized homolog is the of Bacillus subtilis (aBsu"). 
encoded by sigH (formerly spoOH), which is best-known 
for its role in initiating sporulation, an ultimate differen- 
tiation response to starvation, a^si^ directs transcription 
of genes involved in polar septum formation and pro- 
vokes induction of several regulator genes that in turn 
affect expression of signaling pathways or turn on path- 
ways for endospore engulfment (e.g. via the sigma 
factor) [5,6]. assu'^ is also associated with genetic com- 
petence, which enables the uptake of exogenous DNA 
and its assimilation as new genetic information, leading 
to natural genetic transformation. This transient state 
occurs in about 10% of the cells as part of the same 
nutrient depletion response as sporulation. Obsu^ 
increases expression of one of the two peptide phero- 
mones needed for optimal activation of the master regu- 
lator of the competence pathway ComK [7,8]. While 
cJesu'^ is essential for initiating sporulation, its absence 
reduces, but does not abolish transformation (efficiency 
is decreased by ~ 16-fold) [9]. The whole decision-mak- 
ing pathway leading to sporulation or competence is an 
elaborate signal transduction network relying on multi- 
ple partners [7,10]. In addition, aBsu"^ reportedly affects 
expression of about 10% of the genome and was pro- 
posed to be involved in the growth transition to station- 
ary phase [5]. 

The position of aBsu"^ in the tree of a^°-type sigma 
factors is unclear. It exhibits structural characteristics 
similar to ECF sigma factors (group 4), yet phylogenetic 
analyses placed it between groups 3 and 4 [2,4,11]. 
Unexpectedly, database and phylogenetic analyses 
revealed that a'^-like gene products are largely distribu- 
ted in various members of the Firmicutes, sporulating or 
not [11,12]. Interestingly, o -like factors appear to be 
more divergent across non-sporulating bacteria than in 
sporulating bacteria [12]. At the same time, structural 
elements similar to the conserved Gram-positive DNA 
uptake machinery appeared to be encoded in the gen- 
ome in members of the Firmicutes not known for being 
naturally transformable, suggesting that this capacity 
may be more widespread than previously expected 
[12-14]. Two factors, classified in a single large a^- 
family of sigma factors by Morikawa et al. [12], are 
directly involved in transcription of competence genes 
in non-sporulating bacteria: the well-known ComX of 
naturally transformable streptococci [15], and the pro- 
duct of the so-called sigH gene of Staphylococcus aureus, 
a species which has not yet been shown to be trans- 
formable [12]. These observations suggested the link 



between o -like factors and genetic competence in non 
sporulating Firmicutes [12]. 

L. sakei belongs to the microbiota that develops on 
meats under storage, especially during vacuum packing. 
It is largely used as a starter for the manufacture of fer- 
mented sausages in Western Europe and its potential 
use in meat product biopreservation is currently under 
study [16-18]. Survival of L. sakei ranges from one day 
in aerated chemically defined liquid medium, to a few 
months in dry sausages, although little is known about 
the factors determining its stability. The existence in L. 
sakei of sigHi^^ai an apparent sigH^^^ ortholog, led us to 
identify the gene set regulated by a^sa^t and to deter- 
mine whether and how this regulator is implicated in 
competence and stationary phase survival. A strain 
allowing experimental sigHj^^^ induction was con- 
structed, and used in a genome-wide microarray study. 
Genes activated by sigHi^^n overexpression appeared 
mainly involved in genetic competence, although we 
could not obtain evidence for natural transformation. 
This study provides further suggestive evidence that the 
conserved role of the a'^-like sigma factors in non-spor- 
ulating Firmicutes is to activate competence gene 
expression. 

Results and discussion 

Identification of sigH in the genome of L. sakei and other 
lactobacilli 

Automatic annotation of the L. sakei 23 K genome [16] 
identified LSA1677 as a coding sequence (CDS) of a 
putative alternative sigma factor of the superfamily. 
It belongs to COG1595 (E-value of 7e''^), which com- 
prises both ECF-type sigma factors (£. coli RpoE homo- 
logs) and of B. subtilis, and thus reflects the 
reported structural proximity between ECF sigma factors 
and Cbsu" [2,4,11]. The conserved genetic context of the 
L. sakei LSA1677 locus and the B. subtilis sigH locus, 
and more generally the local synteny between several 
members of the Firmicutes (Figure 1), revealed that 
LSA1677 and sigH^su are likely orthologous genes, 
belonging to a widespread family in the Firmicutes. 

Nevertheless, the locus comprising a -like gene may 
have experienced genetic rearrangements across the dif- 
ferent genera and also among species of the same genus 
(Figure 1). Moreover, the a"-like gene location seems to 
be variable in members of the Firmicutes, especially in 
the Lactobacillales (Figure 1). A putative a'^-like gene is 
not found at the same location in Lactobacillus salivar- 
ius as in L. sakei (locus cysS-nusG). Likewise, the loca- 
tion of the unique gene for the ComX factor differs in 
the naturally competent Streptococcus thermophilus 
LMD9 from those of each of the identical comX copies 
in S. pneumoniae R6, in which both copies are adjacent 
to a tRNA gene and ribosomal operons. 
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Figure 1 Genetic context of sigH and a"-like genes in members of the FIrmicutes The figure was generated using IVlicrobes on line 
facilities littp://www.microbesonline.org. Similarly filled arrows represent homologous CDSs. White arrows indicate CDSs without counterpart 
Pseudogene is indicated by a dotted outline. RNA-encoding genes are represented by thin arrows. Two loci are shown for L salivarius, one 
demonstrating the absence of a s/gH counterpart in the same genetic context as 6. subtilis and the other, at a distance of 0.9 Mb, showing the 
sigH homologous gene in its genetic context. Two loci are also shown for 5. pneumoniae, which possesses two identical copies of comX. 
Positions of primers AML50 (upstream) and AML58 (downstream) are indicated by small arrows under the L sakei sigH locus. Species are 
represented by the same strains as listed in Figure 2. 



Although the genetic context of the a -hke locus is 
very well conserved between L. sakei and Lactobacillus 
plantarum, the two a'^-like proteins share only 29% 
amino acid (aa) identity. Indeed, the level of inter-spe- 
cies aa identity of a'^-like gene products across the 
genus Lactobacillus is low (e.g., < 20% between L. 
plantarum WCFSl and L. jensenii 208-1 to 67% 
between L. helveticus DPC4571 and L. crispatus 
MVIAUS). The LSA1677 gene product shares weak aa 
identity with the a" factors of B. subtilis (24%) and S. 
aureus (21%), as well as 22% aa identity with ComX of 
S. pneumoniae (see Additional file 1: Alignment of 
four a'^-group sigma factors). Due to the high 
sequence divergence between sigma factors, a robust 
phylogeny is difficult to achieve. Tentative clustering 
of a'^-like sigma factors (Figure 2), also including spor- 
ulation and known ECF sigma factors of B. subtilis, 



separates agsu from the other sigma factors in that 
species and argues for the existence of a a'^-type 
family in Firmicutes [12]. a^-like factors appear to 
form groups mostly congruent with the genus phylo- 
geny, irrespective of the location of the relevant gene 
in the genomes (Figure 2). The a"-like sigma factors 
of lactobacilli added a fourth group to the three pre- 
viously reported groups (whose type factors are Gbsu". 
a'^-like of staphylococci and ComX of streptococci) 
[12]. LSA1677 exhibits the characteristics of a a'^-like 
factor encoding gene and, in view of the conserved 
genetic context between L. sakei and B. subtilis, was 
called sigH. Note that the name sigX has been chosen 
for recently annotated genomes of Lactobacillales. 
Although the name SigX is more appropriate than 
ComX for a sigma factor, it adds confusion with the 
existing SigX sigma factor of B. subtilis, which is not 
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Figure 2 Clustering of selected (i™-superfamily of sigma factors. The unrooted tree resulted from a multiple alignment over the whole aa 
sequence length of a^-like factors and known sigma factors from group 3 (sporulation factors of S. subtilis) and group 4 (ECF factors from 6. 
subtilis and Gram-negative bacteria). The multiple alignment was generated using clustalX [19], by introducing first the shortest sequences to 
ensure a correct alignment of the conserved regions. The tree was drawn with NJplot http://pbil.univ-lyon1.fr/software/njplothtml. Bootstrap 
values (number of seeds: 1000, number of trials: 100) are indicated for the upper branches. Evolutionary distance is represented by branch length 
(scale at the bottom). Groups of a"-like factors were numbered as previously reported [12] and a fourth group (IV) was added by our analysis. 
Bsu, Bacillus subtilis 168; EC, £ coll K-12 substr. MG1655; Pae, Pseudomonas aeruginosa PA01; Ef, Enteroccocus faecalis V583; Lla, Lactococcus lactis 
11403; Lmo, Listeria monocytogenes EGD-e; Genus Clostridium: CBO, C botulinum A ATCC3502; CP, C. difficile 630. Genus Lactobacillus: Lba, L 
acidophilus NCFM; Lsei, L casei ATCC334; Lgas, L. gasseri ATCC 33323; Lp, L. plantarum WCFS1; Lsa, L. sakei 23 K, Lsl, L salivarius UCC1 18; Lac, L 
acidophilus NCFM. Genus Staphylococcus: Sau, S. aureus N315; Sea, 5. carnosus TM300; SE, S. epidermidis ATCC 12228. Genus Streptococcus: Spn, S. 
pneumoniae R6; Spy, S. pyogenes ATCC 10782; Sth, S. thermophilus LMD-9. Names of gene products or locus tags are indicated, a^-like sigma 
factors which belong to sporulating bacteria are indicated with an asterisk; those encoded by a gene not located at a similar locus to s/gHesu are 
underlined (dashed line for the particular case of S. pneumoniae, see Figure 1). The best studied a^-like sigma factor for each group is in bold 
type. 



the equivalent of a . This certainly calls for a unified 
nomenclature of sigma factors in Firmicutes. 

Conservation of sigH genes in the L sakei species 

We asked whether sigH genes were conserved among L. 
sakei isolates exhibiting a broad intraspecies diversity 
[50]. Based on the presence or absence of markers of 
the flexible gene pool, L. sakei isolates from various 
sources were previously classified into distinct genotypic 
clusters, possibly affiliated with two prevailing sub- 



species [20]. The 5' and 3' ends of the sigH gene were 
used as targets for PGR amplification of 17 isolates 
belonging to 9 of the 10 reported clusters of the species 
[20]. A unique fragment of the same size as that of 
strain 23 K was amplified for all samples, indicating that 
the gene is likely present in all tested isolates. Further 
analysis of the locus was undertaken for 7 of these 
strains distributed in 5 clusters. Amplification obtained 
with primers designed on the basis of the L. sakei 23 K 
genome outside of sigH suggested that the genetic 
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context is conserved in all these strains (see position of 
primers AML50 and AML58 in Figure 1). Polymorphism 
analysis of the sigH sequences brought additional infor- 
mation. As shown in Figure 3, 29 polymorphic sites 
were identified in the sigH CDS, of which only 9 were 
involved in 7 aa changes, mostly conservative. Thus, 
SigH function and coding gene location appear to be 
conserved in the L. sakei species. 

L. lactis subspecies lactis and cremoris exhibit two 
comX allelic types whose nucleotide divergence is at 
most 27.5% [21]. In contrast, sigH divergence (4.5% 
maximum divergence) was incongruent with the pre- 
viously reported genotypic classification of L. sakei 
strains [20], and its two proposed subspecies (Figure 3). 
This discrepancy may be explained either by a particular 
evolutionary history of that gene in L. sakei or by the 
possibility that the classification based on the flexible 
gene pool does not reflect the phylogenetic relationships 
between strains which remain to be established. 

High nucleotide divergence between species, con- 
trasted with generally higher conservation within spe- 
cies, was also observed for sigH loci in the genus 
Staphylococcus [22]. The reason for such high inter-spe- 
cies polymorphism is unknown. However, rapid evolu- 
tion after species divergence rather than lateral gene 
transfer may be responsible, as the phylogeny of sigH 
genes was reported to be concordant with species phylo- 
geny in staphylococci [22]. 

As reported in this paper, functional studies were further 
conducted on RV2002, a derivative of L. sakei strain 23 K, 
for which genome data is available, and in which the endo- 
genous P -galactosidase encoding gene is inactivated, thus 
enabling the use of a lacZ reporter gene [23]. 



PCR (qPCR) following reverse transcription to determine 
if sigHi^sa expression in L. sakei is also temporally regu- 
lated. L. sakei was cultivated in chemically defined med- 
ium (MCD) at 30°C and total RNA was extracted from 
cells 2 h after inoculation and every hour from 4 to 8 h. In 
these experiments, transition to stationary phase was 
observed between 5.5 and 7 h. Three genes, Idh, gyrA and 
sigA, were initially evaluated as candidate internal stan- 
dards for qPCR, based on previously used standards in 
Oenococcus oeni [25]. We selected Idh, which showed the 
least variation of mRNA levels during growth (Figure 4). 
sigHi^sa mRNA levels were then quantified relative to the 
early-exponential condition (2 h) chosen to calibrate the 
measurements, and by normalizing with Idh mRNA. 
Results showed a slight increase (1.7 ± 0.3) of sigHi^sa tran- 
scripts around the transition to stationary phase (Figure 
4). This transcription pattern is close to that reported for 
B. subtilis, for which sigH^^u transcription reached a 3-fold 
increase peak 40 min before transition to stationary phase 
in sporulation medium [24]. Possibly, the observed level of 
sigHi^s^ induction could be greater in other media and 
growth conditions. sigH-g^^^ repression during exponential 
growth phase relies on the transcriptional repressor AbrB, 
a major transition-state regulator in B. subtilis [24]. As no 
homolog of AbrB could be identified in L. sakei, we sus- 
pect that other regulatory circuit may be involved in con- 
trolling sigHi^ssL- Interestingly, S. aureus sigH^^^ 
transcription reportedly decreases 10-fold from early- 
exponential to stationary phase [26]. 

Overexpression of 

The sigH]^^^ gene was overexpressed as a means to 
reveal genes that it specifically regulates. sigHj^^a was 



Temporal transcription of sigH 

In B. subtilis, sz^/Zbsu transcription increases from mid- 
exponential to stationary phase [24] . We used quantitative 
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Figure 4 Temporal transcription of sigH. Growth of R\/2002 has 
been monitored by ODgoo (right axis). Time is indicated in hours 
relative to the approximate transition to stationary phase (T). mRNAs 
levels of Idh (grey blocks) or sigH (white blocks) were measured by 
qPCR and expressed as fold change relative to an early-exponential 
calibrator sample (left axis). For sigH, results have been further 
normalized by idh mRNA levels and expressed as sigH/idh ratio. 
Error bars represent standard deviation. A fold change of 1 indicates 
a constant level of transcripts. 
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placed under the control of the copper-inducible L. 
sakei promoter PatkY) present on plasmid pRV613 [27], 
and the resultant plasmid was introduced into RV2002 
wild-type (WT) strain. The resulting strain, designated 
sigH{hyY, thus has an additional expression-controlled 
copy of sigH and was compared to the equivalent WT 
strain harboring the pRV613 plasmid, in which PatkY 
controls lacZ (see additional file 2: Genotype of L. sakei 
strains affected in sigH). We anticipated that compe- 
tence genes, found in the L. sakei genome and likely 
coding for a DNA uptake machinery [28], might be tar- 
get genes for transcription by a"-directed RNA poly- 
merase (see additional file 3: Competence DNA uptake 
machinery of B. subtilis and comparison with L. sakei). 
We therefore used qPCR to monitor sigHi^^.^ and comEA 
transcripts produced by cells at 1 h (the time needed to 
obtain ~15-fold induction of Patky! [27]) and 4 h (corre- 
sponding roughly to the transition to stationary phase) 
after induction in mid-exponential phase, or without 
induction. Addition of CUSO4 to the strain harboring 
the control plasmid had no detectable effect on the 
amount of sigHj^sa E^nd comEA transcripts (Table 1). In 
contrast, induction of the Patky-controlled copy of 
sigHi^sa led to a ~40-fold effective increase of sigH tran- 
scripts after 1 h, and ~ 200-fold after 4 h. comEA tran- 
script levels were highly increased (over 300-fold), but 
only when sigHi^^a was 40 fold over-expressed (a 20-fold 
increase of sigHj^nn mRNA did not alter comEA expres- 
sion. Table 1). The need for high sigHi^sa overexpression 
may indicate the need to overcome posttranscriptional 
controls to produce enough active CTLsa"^- This proposal 
is supported by observations in B. subtilis, where a^su 
was shown to be subjected to multiple controls [5,29], 
and in the genus Streptococcus, where high levels of 
ComX overexpression were required to artificially 
induce competence [30], likely due to the negative con- 
trol of ComX stability by a Clp protease complex 
[30,31]. 

comEA transcription was not increased at the onset of 
stationary phase in the WT nor in the induced sigHihy)" 



strain, suggesting that the competence genes are not 
naturally induced under laboratory conditions. Activa- 
tion of comEA tended to diminish after a four hour- 
induction despite high levels of sigHi^sa transcripts, pos- 
sibly indicative of another regulatory loop on comEA or 
post-transcriptional regulation of sigHi^sa- This transcrip- 
tion pattern was similar for comGA exhibiting a 280- 
fold increase in transcript amounts one hour after 
sigHi^s^ induction in sigH{\\y)" followed by a 3-fold 
decrease between one and four hours. These results 
show that in L. sakei, conditions of Glss overexpression 
lead to activation of comEA and comGA. Nevertheless, 
other factors likely modulate com gene expression, as 
suggested from the drop of expression late in growth. 

Global transcriptome profiling of the overexpression 
strain suggests the existence of a genetic competence 
state 

To identify the genes of L. sakei regulated by aLsa"^; the 
experimental system described above was used in a full- 
genome comparative transcriptome analysis of s/gW(hy)* 
and sigH{wt)" after one hour induction with 30 [iM 
CUSO4. Quantification and statistical analysis of the 
microarray data (see Methods for parameters) led to 
relatively few differentially expressed candidate genes. 
The overexpressed sigH gene in sigHihy)'' was 11 ± 3 
times induced compared to the WT strain in this micro- 
array experiment; qPCR-based quantification of the 
same RNA samples showed a 149 ± 42-fold greater 
expression relative to the WT strain, confirming the 
successful overexpression of sigHi^^n- Differences in fold 
ratios between microarray-profiling and qPCR analysis 
are not unusual but were high in our experiment; they 
might reflect a less efficient detection on microarray or 
an overestimation by qPCR especially when genes are 
weakly expressed in one of the conditions, which 
seemed to be the case for the com genes. 

Based on statistical tests (P value < 0.05), our microar- 
ray analysis initially identified some 25 candidate genes 
whose expression was likely affected by sigHj^^^ 



Table 1 Relative expression ratio^ of sigH and comEA with or without overexpression of sigH 



Sample 


sigH{\Nt)* i 


s/gH(hy)* ni 


s/gH(hy)* 1 


Calibrator 


s/gH(wt)* ni 


s/gH(wt)* nl 


s/gH(wt)* i 


Measured effect 


control of the Inducer effect in wt 
strain 


presence of the additional copy of 
sigH 


cumulative effect of Induced additional 
copy 


Time (h) 


1 4 


1 4 


1 4 


sigH 


1 1 


7 24 


40 200 


comEA 


1 1 


2 3 


370 80 



Expressed as fold change of transcripts amounts of each gene in each given sample relative to the indicated calibrator and normalized with Idh. Results are the 
mean of two independent experiments. The level of Idh transcripts was stable, irrespective of the copy number or induction status of sigH (e.g. mean fold 
change across all induced samples relative to non induced samples: 0.9 ± 0.2). Note that sigH is present at one (chromosomal) copy in sigH(wt)* and at two 
copies (one additional copper-controlled copy on a plasmid) in s;gH{hy}*; the transcription of both is measured simultaneously, ni and i refer to 'not induced' and 
'induced', respectively. 
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overexpression; behavior of several genes was confirmed 
by qPCR (Table 2). The known genes can be grouped 
into two main functional categories: competence (DNA 
uptake) and DNA metabolism. All the late competence 
(com) operons encoding structural elements of the DNA 
uptake machinery were highly activated by sigHi^^a over- 
expression. In contrast, transcription of ssb, regulated as 
a late competence gene in B. subtilis [32], was nearly 
constant or only very weakly induced. Other genes 
involved in DNA metabolism, and known to be induced 
during the competence state in other species, i.e., 
recombination genes recA and dprA, both of which are 
involved in natural bacterial transformation in different 
species [33], gave a contrasted picture when their tran- 
scription was specifically measured by qPCR. Whereas 
recA was little activated, expression of dprA was highly 
induced in the sigH{hyY context (Table 2). 

Two genes coding for hypothetical proteins, LSA0409 
and LSA0779, were down-regulated in the sigHi^sa over- 
expression strain. As sigma factors are usually positive 
regulators, we consider it likely that down-regulation of 
these genes is an indirect effect of sigHi^^;.^ overexpres- 
sion, e.g., this effect could correspond to a -mediated 
activation of an unidentified repressor. The sole tran- 
scriptional regulator (LSA0421) listed as a"-activated in 
Table 2 is probably not responsible for this effect, since 
MerR-type regulators reportedly act as activators [34]. 
No clue for the function of LSA0421 could be drawn 
from the genetic context or multiple alignments with 
other MerR-type proteins from E. coli and B. subtilis. 

Overall, the qPCR analysis validates our statistical ana- 
lysis of the microarray data based on the common var- 
iance model associated with the correction of 
Bonferroni (Table 2). Indeed, although calculated 
expression ratios were very similar for comFA and ssb 
(around 1.5), the former, which had an associated P 
value < 0.05 with the Bonferroni correction, was con- 
firmed as overexpressed in the qPCR analysis, whereas 
the latter (which passed the other applied statistical test) 
was found to be almost unaffected by a" in the qPCR 
analysis (Table 2). Therefore, we expect that all genes 
with a better score than comFA in the microarray anlay- 
sis (i.e. P value Bonferroni < 1.54 E-02) are good candi- 
dates for belonging to the aLsa" regulon. Altogether, 
results of this study thus identify 25 genes as belonging 
to the aLsa"^ regulon. Some genes (e.g., dprA), while 
truly activated by CLsa'^. may not be detected in this 
microarray experiment, indicating the need for further 
studies to define the full regulon. 

Transcriptional reprogramming caused by sigH]^^^ 
overexpression is consistent with the existence of a 
competent state in L. sakei, supported by the observed 
up-regulation of com genes involved in pseudopilus 
morphogenesis and DNA translocation as well as of 



dprA (which shares 47% aa identity with the 5. pneumo- 
niae dprA gene product), ssb and recA appear little or 
not activated one hour after sigHj^sn induced overexpres- 
sion, whereas their level of induction during the compe- 
tence state in S. pneumoniae and B. subtilis reportedly 
varies from 5 to over ten-fold [32,35]. These genes 
might be transiently regulated (in a narrower window 
than com operons and dprA), regulated by other factors, 
or their up-regulation may not be required in L. sakei. 
Indeed both genes participate in the bacterial vegetative 
life cycle and are expected to be transcribed at a signifi- 
cant basal level when cells are not in the competence 
state [36]. Interestingly, L. sakei possesses a unique ssb 
gene (ssM-type), whereas B. subtilis and S. pneumoniae 
have paralog genes [36,37]. The need for a transforma- 
tion-dedicated SSB protein has been discussed [37]. 
Although known natural transformation is frequently 
associated with multiple ssb in Firmicutes [37], an addi- 
tional competence-induced SSB may be a facilitator 
rather than an essential contributor to the transforma- 
tion process, since transformation frequency is only 
reduced by two- to ten-fold when ssbB is inactivated in 
S. pneumoniae or in B. subtilis [36]. 

Is L sakei capable of natural genetic transformation? 

As the a -activated transcriptome of L. sakei was indi- 
cative of a competence state, we looked for genetic 
transformation in this species. The first strategy involved 
the overproducer strain sigH{\\y)". As the strain harbors 
a resident plasmid conferring erythromycin resistance, a 
compatible plasmid conferring chloramphenicol resis- 
tance (pVI1056) was used for transformation tests. After 
cultivation under inducing conditions (i.e., addition of 
30 [iM CUSO4), the strain was mixed with 100 ng of 
pVI1056 and plated on selective medium. Experiments 
were performed under various conditions: i) glucose 
concentration at 1% or 0.1%, ii) growth in microaerobio- 
sis or aeration, and induction at early, middle or late 
exponential phase iii) addition of MgCl2 (80 mM) dur- 
ing contact between cells and DNA, after middle phase 
induction in microaerobiosis or aeration; in addition, 
chromosomal L. sakei DNA (1 |ig) was also used as exo- 
genous DNA. None of the tested conditions resulted in 
DNA transformation. 

Development of natural transformation may be strain- 
dependent [30,38,39]. We therefore used a second strat- 
egy (independent of sigH overexpression) to test differ- 
ent L. sakei isolates for competence, using a protocol 
where DNA and strains are deposited on solid medium. 
In addition to 23 K, four strains (64 K, 332 F, 160 K 
and LTH675) were chosen based on their different gen- 
otypes and genome sizes, and known capacities to be 
transformed by electroporation [20,58]; Chaillou and 
Anba, personal communication]. Two replicative 
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Table 2 Genome-wide transcriptome profiling of SigHLsa overexpression in L.sakei 23 K 



Functional Gene 
category and CDS Name 



Product 



Rvalue 
(Bonferroni) 
common 
variance 
model 



Rvalue (FDR) 
varmixt model 



Expression 
s/gH(hy)*/ 



microarray 



ratios 
s/gH(wt)* 



qRCR 



Competence 



LSA0492 


comFA 


DNA 


uptake machinery § 


1 .54E-02 


> threshold 


1.5 ± 


0.4 


LSA0493 


comFC 


DNA 


uptake machinery 


0 


3.56E-03 


2.2 ± 


0.2 


LSA1069 


comEC 


DNA 


uptake machinery 


9.52E-10 


1.31 E-02 


1.9 ± 


0.2 


LSA1071 


comEA 


DNA 


uptake machinery 


0 


7.23E-03 


2.5 ± 


0.3 


LSA1301 


comGF 


DNA 


uptake machinery 


0 


2.71 E-04 


3 ± 


2 


LSA1302 


comGE 


DNA 


uptake machinery 


0 


1 .44E-06 


3.7 ± 


0.5 


LSA1303 


comGD 


DNA 


uptake machinery 


0 


2.21 E-04 


2.8 ± 


0.3 


LSA1304 


comGC 


DNA 


uptake machinery 


0 


5.62E-12 


7 ± 


2 


LSA1305 


comGB 


DNA 


uptake machinery 


1.02E-10 


3.57E-02 


2.0 ± 


0.3 


LSA1306 


comGA 


DNA 


uptake machinery 


3.17E-09 


7.25E-03 


1.9 ± 


0.2 


LSA1771 


comC 


DNA 


uptake machinery 


0 


410E-06 


3.2 ± 


0.2 



DNA metabolism: replication, repair, recombination, RM 



261 



115 



421 ± 
104 



199 



LSA0008 ssb 


Single-stranded DNA binding protein 


> threshold 


3.88E-02 


1.4 ± 


0.1 


1 .2 ± 0.3 


LSA0146 


Putative DNA methyltransferase (apparently stand- 
alone) 


1.55E-04 


> threshold 


1.6 ± 


0.4 




ISA 1299 


Putative DNA methyltransferase (apparently stand- 
alone) 


248E-08 


> threshold 


1.9 ± 


0.4 




ISA 1338 exoA 


Exodeoxyribonuclease II 


1 .36E-07 


> threshold 


1.8 ± 


0.3 




Purines, pyrimidlnes, nucleosides and nucleotides 












LSA0533 iunh2 


Inosine-uridine preferring nucleoside hydrolase 


1.14E-05 


> threshold 


1.7 ± 


0.4 




Energy metabolism 














LSA1298 ack2 


Acetate kinase 


4.27E-09 


> threshold 


1.9 ± 


0.4 




Translation 














LSA0009 rpsR 


Ribosomal protein 


1 .67E-02 


> threshold 


1.5 ± 


0.4 




Regulatory function 














LSA0421 


Putative transcriptional regulator, MerR family 


0 


3.56E-03 


2.5 ± 


0.5 




Hypothetical protein 














LSA0040 


Hypothetical protein, conserved in some lactobacilli 


0 


3.56E-03 


2.5 ± 


0.5 




LSA0409 


Hypothetical integral membrane protein 


3.02E-05 


7.25 E-03 


0.61 ± 


0.01 




LSA0536 


Hypothetical protein with putative NAD-binding 
domain, NmrA structural superfamily 


6.28E-06 


3.32E-02 


1.6 ± 


0.4 




LSA0779 


Hypothetical protein, peptidase S66 superfamily 


477E-05 


> threshold 


0.6 ± 


0.1 




LSA0991 


Hypothetical protein with putative NAD-binding 
domain, NmrA structural superfamily 


1 .02E-04 


> threshold 


1.6 ± 


0.2 




LSA1475 


Hypothetical protein, conserved in bacteria 


1.62-12 


> threshold 


2.1 ± 


0.5 




CDS £ Gene 


Product 










qPCR 


Name 














LSA0487 recA 


DNA recombinase A 










2.7 ± 0.7 


LSA0992 dprA 


DNA protecting protein, involved in DNA 
transformation 










2163 ± 
1242 



$ Expression ratios represent the fold change in amounts of transcripts in the strain overexpressing SigH relative to the WT control strain. For the microarray 
experiment they were calculated from log2ratio; for the qPCR they were calculated by the 2"'^^'"^ method described in Methods. Genes underexpressed in the 
context of SigH overexpression have a ratio < 1. Standard deviation is indicated {weak accuracy for qPCR experiments may be due to Ct at the detection limit for 
basal level). 

§ see additional file 3: Competence DNA uptake machinery of B. subtilis and comparison with L. sakei. 

£ not found statistically differentially expressed in the microarray transcriptome experiment, checked by qPCR. 
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plasmids and chromosomal L. sakei DNA were used. In 
spite of varying media (MRS or MCD) and incubation 
temperatures (4°C, 30°C or 37°C), no colonies were 
recovered on selective medium. 

Among the Lactobacillales, natural genetic transfor- 
mation has been reported for many species of the genus 
Streptococcus [40] and has been suspected for one Lac- 
tobacillus [41]. In recent years, natural transformation 
has been demonstrated in several Gram-positive or 
Gram-negative species, previously unsuspected to 
develop genetic competence [42,43]. Overproduction of 
the activator protein has proven to be an efficient way 
to trigger genetic transformation in various bacteria, e. 
g., TfoX in Vibrio choleras [42] or ComK in Bacillus 
species [14,44]. However, artificially raising transcription 
of the ComX master regulator gene initially failed to 
induce efficient genetic transformation for S. thermophi- 
lus strain LMD-9 [30], which was very recently shown 
to be efficiently naturally transformable [37]. In the pre- 
sent and previous studies, a failure to achieve a compe- 
tent state in bacteria (either spontaneously or triggered 
by artificial overexpression of a master activator) may be 
due to the use of inappropriate growth conditions, 
which might not allow the detection by the cells of a 
needed specific triggering factor [38,42] or the full acti- 
vation of multiple converging regulatory pathways [30]. 
As such, in the case of L. lactis [21], S. pyogenes [45], S. 
aureus [12], or L. sakei (this paper), only the activation 
of several competence genes, but not genetic transfor- 
mation, could be obtained after ectopic expression of 
the activating sigma factor. Our results suggest that 
some of the genes induced in other naturally competent 
Firmicutes are not activated by the sole sigHi^sa overex- 
pression in L. sakei. In the case of 5. pyogenes, the iden- 
tification of a novel pheromone in related species of 
Streptococcus might pave the way for deciphering a nat- 
ural genetic transformation system in this bacterium 
[46]. Whether competence gene activation by ComX/o 
is linked to the capacity of being transformable in these 
species, and under which conditions, remains to be 
determined. 

Effect of sigH on L. sakei survival 

No indication of another large adaptive response trig- 
gered by aLsa'^ could be deduced from the few other 
up-regulated genes distributed in different functional 
categories. We also searched for phenotypic effects 
linked to a putative role of a" on survival in stationary 
phase or after DNA damage. For that purpose, we con- 
structed a sigH(nul) null mutant (see Methods) and 
compared the effect of overexpression or absence of 
t^Lsa'^ relative to WT strains on growth and stationary 
phase survival in MCD medium under aerobiosis, 
microaerobiosis or anaerobiosis, as well as on UV 



resistance. No changes in any of the above tests could 
be attributed to a expression levels under the condi- 
tions tested (data not shown). Interestingly, all the 
strains revealed UV resistance, since the fraction of each 
population killed by 254 nm irradiation was in the range 
of 0-5% at 60 J.m 60-70% at 80 J.m ^ 95-98% at 100 J. 
m'^ and 99.5-99.9% at 120 J.m'^. This is to be compared 
to the reported 100% killing of Lactobacillus brevis 
exposed to 254 nm UV light at 70 J.m'^ [47]. Competi- 
tion experiments in mixed cultures revealed no imbal- 
ance in growth or survival between the a 
overproducing or a deficient and WT strains in MCD 
medium (Figure 5). As MCD medium may not represent 
a usual environment for the bacterium, a meat-derived 
medium was tested for comparison of sigH{nul) and 
WT strains. L. sakei showed prolonged stationary phase 
survival in meat juice, where about one percent of the 
population was still alive after one month at 30°C (Fig- 
ure 6). Inactivation of sigH brought no striking change 
to the phenotype. 

Conclusions 

This study gives further insight into the function of a^- 
family sigma factors from Firmicutes, whether they 
belong to sporulating or non-sporulating bacteria. Two 
main models emerge: ComX of non-sporulating strep- 
toccoci is the master activator of the genetic compe- 
tence pathway, participating directly in the transcription 
of late com genes [15]. a of B. subtilis activates a com- 
plex response leading to spore formation as an ultimate 
outcome and to the development of genetic competence 
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days 



Figure 5 Effect of overexpression or deletion of sigH on 
viability of L sakei in mixed cultures with WT strain Each pair 
of mutant and WT strains has been mixed after separate growth 
until an OD500 of 0.3, in MCD medium at 30°C in microaerobiosis. 
Enumeration on appropriate agar plates allowed to distinguish WT 
from mutant strains. s/gH(nul) mutant (black triangles) was mixed 
with V\^ strain 23 K (empty triangles). sigH{hy)* overexpression 
mutant (black circles) was mixed with sigH(wt)* strain (empty 
circles), and 30 |jM CUSO4 was added to the culture. Curves are the 
mean of two independent experiments. 
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Figure 6 Long-term viability of L sakei in meat juice at 30°C 

Curves are the mean of three independent experiments; error bars 
represent standard deviation (logarithmic scale). 



during a transition period. Unlike ComX, a^su does not 
directly activate genes encoding the DNA uptake 
machinery, but participates as an intermediate in the 
upstream signaling pathway controlling the master regu- 
lator of competence ComK [5,48]. sigH genes from the 
non-sporulating L. sakei and S. aureus species are orga- 
nized similarly to the sigH locus of the sporulating bac- 
terium B. subtilis. However, unlike B. subtilis, they act 
like streptococcal ComX by activating late com genes 
[[12]; this paper]. We speculate that this function may 
be conserved in the order Lactobacillales, irrespective of 
the exact location of the so-called ComX or a encod- 
ing gene. The regulon of CLsa as deduced by assessing 
the effects of CLsa overexpression was rather small. It 
should be mentioned that the genome size of the model 
strain used was 136 kb less than the average size within 
the species [20] and that our strategy mainly identified 
genes that were strongly affected by aLsa," indepen- 
dently of possible other, undetermined, environmental 
signals. A large number of reported regulatory effects of 
CTbsu"^ are actually mediated in conjunction with other 
transcriptional regulators, especially SpoOA and AbrB 
[5]. L. sakei and more generally Lactobacillales do 
apparently not possess orthologs of these regulatory 
proteins, neither do they possess a ComK homologue. 
Deciphering all the functions of the conserved sigma 
factor in other groups of Firmicutes, sporulating or not, 
and equipped with different combinations of these 
known global regulators will probably help to clarify a 
evolution in this group of bacteria. 

Methods 

Media and growth conditions 

L. sakei was grown at 30°C in MRS medium [49] or in 
the chemically defined medium MCD [50], both con- 
taining 1% glucose. A two-step preculture was used to 
assure reproducibility of experiments. First, 5 ml MRS 
was inoculated with one freshly isolated colony and 
incubated for about 8 h without agitation. After 



centrifugation, cells were resuspended in MCD at an 
ODgoo of 1 and 10 to 20 [il of the suspension was used 
to inoculate 40 ml of fresh MCD. This second precul- 
ture was incubated without agitation for about 15 h so 
as to collect the cells in exponential growth phase. This 
preculture was then concentrated to an ODgoo of 10 in 
fresh MCD, and used to inoculate the test culture to 
give an initial ODgoo of 0.1 to 0.15. Unless otherwise 
indicated, growth conditions under microaerobiosis 
were used. Different aeration conditions were obtained 
by varying the agitation parameter and volume of cul- 
tures. Typically microaerobiosis corresponded to non- 
agitated cultures in closed centrifugation tubes, aerobio- 
sis to 140 rpm of 1/10* full ehrlenmeyer flasks, and 
anaerobiosis to N2-saturated cultures in sealed bottles 
[51]. Growth was followed by ODgoo measured in a 
Secomah spectrophotometer. As 30 [^M CUSO4 may be 
added to the culture, we monitored its global effect on 
L. sakei growth. In static or anaerobic growth condi- 
tions, 30 |iM CUSO4 had no effect on growth. In aera- 
tion conditions, 30 \xhA CUSO4 had a slight effect on 
growth (2-10% lower ODgoo at the end of growth), and 
slightly extended viability. 

Meat juice was obtained from beef meat homogenized 
with half volume of sterile water in a Stomacher for 2 
cycles of 3 min each. The supernatant obtained after 
centrifugation (10,000^ for 15 min) was filter sterilized 
and stocked at -20°C (M.-C. Champomier Verges, 
unpublished). Escherichia coli (DH5aF' or TGI) was cul- 
tured aerobically in LB at 37°C. Selective pressure for 
plasmids was maintained in E. coli with ampicillin 100 
mg.r'^, and in L. sakei, with erythromycin 5 mg.l'^. 

DNA techniques 

Standard procedures were used for DNA manipulation. 
Classical PCR reactions were performed with Taq poly- 
merase (Fermentas) or Pfu polymerase (Promega) for 
cloning purpose, and run in MJ research PTC-200 ther- 
mocycler. Extraction of plasmids and chromosomal 
DNA as well as electroporation of L. sakei and L. casei 
BL23 was carried out as described [52]. Primers are 
listed in additional file 4. 

Diversity of sigH in L sakei 

L. sakei strains (18, 21, 23 K, 64, 112, 160 K, 300, 332, 
JG3, MF2091, MF2092, ATCC15521, CIP105422, SF771, 
LTH677, LTH2070) were from our collection or differ- 
ent sources as described [20]. PCR amplification of the 
sigH locus was carried out with two pairs of primers 
(AML31/AML32 and AML50/AML58). Sequence of the 
561 nt fragment corresponding to entire CDS and the 
77 nucleotides of the upstream intergenic region was 
performed on PCR-amplified genomic DNA using each 
of the four primers. Pairwise distances were calculated 
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by MEGA 4 [53] using a Kimura 2-parameter substitu- 
tion model. 

Construction of sigH mutant and sigH expression strains 

SigH production and sigH mutant strains were con- 
structed from RV2002, a derivative of L. sakei 23 K that 
had undergone a deletion of the lacLM gene encoding 
P-galactosidase [23]. Their construction used plasmids 
pRV610 and pRV613 [27] which contain two replication 
origins, one functional in E. coli (pBluescript) and one 
for Gram-positive bacteria (pRV500). 

The L. sakei overproducer strain sigHiYiy)" was 
obtained by introducing plasmid pRV619 into RV2002. 
pRV619 was constructed from pRV613 which bears the 
PatkY copper-inducible promoter cassette of L. sakei 
fused to the E. coli lacZ reporter gene [27]. lacZ was 
replaced by sigHi^^sL in pRV619 as follows. The sigHi^^a 
coding region was PCR-amplified from L. sakei strain 23 
K chromosomal DNA with primers AML31 and AML32 
and the BamHI/Xbal fragment was cloned into pRV613 
digested by the same enzymes, using Lactobacillus casei 
BL23 as a host, since neither L. sakei nor E. coli were 
successful for direct cloning. The extracted plasmid was 
then successfully introduced into E. coli from which we 
could obtain the needed high-quality plasmid prepara- 
tion (Qiagen Plasmid Kits) to electroporate L. sakei 
RV2002. 

The L. sakei sigH null mutant (RV7003 designated 
sigH{n\x\) in the text) was obtained via a double cross- 
over homologous recombination with the pRV622 inte- 
grative plasmid. To inactivate the sigH gene we deleted 
its putative promoter and the first 34 codons while 
introducing an in-frame stop codon at the endpoint of 
the deletion (see additional file 2: Genotypes of L. sakei 
strains affected in sigH). The upstream and downstream 
fragments were generated by PGR using respectively 
AML51/AML52 and AML53/AML54 primer pairs, 
thereby introducing an EcoRI site in sigH. Each ampli- 
con was digested with EcoRI, followed by DNA ligation 
and digestion with PstI and Xhol. The resulting 1.1 kb 
fragment was then reamplified by the distal primers 
AML51 and AML54 and cloned by blunt-end ligation 
after treatment with T4 polymerase, into the pRV610 
cloning vector [27] cut by Smal. As above, L. casei BL23 
was used as a host for cloning, giving plasmid pRV621. 
This plasmid was then successfully introduced into E. 
coli and an intra-molecular deletion of the Gram + 
replication cassette was generated between unique 
restriction sites EcoRV and Kpnl repaired by T4 poly- 
merase, giving pRV622 which replicates in E. coli. Gene 
replacement in L. sakei was carried out as described 
[23], with two successive single crossovers, the first one 
leading to chromosomal integration of the plasmid 
(maintained by erythromycin selection), and the second 



one allowing plasmid excision, monitored by loss of ery- 
thromycin resistance. The mutant chromosomal struc- 
ture was checked by PGR. Correct sequence of the 
inserts was checked for pRV619 and pRV622. 

Induction of Patkv promoter utilization and monitoring 
using ^-galactosldase activity 

The copper-inducible PatkY promoter was used as 
described [27] for sigH overexpression. For this purpose, 
GUSO4 was added to a final concentration of 30 i^M 
when cultures reached an ODgoo of about 0.4. Induction 
of the PatkY promoter was controlled with the sigHiwt)" 
strain, harboring a PatkY-directed lacZ reporter gene. 
Sampling was done one hour after induction and P- 
galactosidase activity was measured according to [23] 
using ONPG (o-nitrophenyl-P-D-galactopyranoside) as a 
substrate. Activities expressed as Miller units relative to 
ODgoo of the culture [23] were observed to be between 
10 and 25 after induction, whereas the non induced 
standard was around 0.5. 

Extraction of total RNA 

L. sakei strains were cultivated at 30°C in MGD under 
microaerobiosis following the standardized procedure 
described in the upper section, in the presence of ery- 
thromycin for plasmid-containing strains. Cultures of L. 
sakei were distributed in as many centrifugation tubes 
as scheduled collecting points and were incubated at 30° 
G without agitation. Gell pellets were collected by brief 
centrifugation at the time indicated, freezed in liquid 
nitrogen and preserved at -80°G. Total RNA was 
extracted with TRIzol reagent (Invitrogen) as previously 
described [54] . Integrity of RNA was checked by Bioana- 
lyzer 2100 (Agilent). RIN values were above 9. 

Whole-genome microarray analysis 

The L. sakei microarray http://migale.jouy.inra.fr/sakei/? 
q=supplement comprises all the identified coding genes 
of strain 23 K represented by 70 nt long oligonucleo- 
tides synthesized by Operon Biotechnologies Inc. The 
manufacture of DNA chips as well as labelling, hybridi- 
zation and image analysis were performed at the Bio- 
chips platform of Toulouse-Genopole http://biopuce. 
insa-toulouse.fr/Maquette/en/. Each oligonucleotide was 
spotted in triplicate on UltraGaps coated slides (Corn- 
ing® Life Sciences). Total RNA (5 |ig) was reverse tran- 
scribed and labeled with either Cy5 dCTP or Cy3 dCTP 
(Amersham Biosciences) using the ChipShot™ Direct 
Labeling System (Promega). Labelled cDNA (50 pmol of 
Gy3 and 50 pmol of Cy5) was included in a dye-switch 
hybridization protocol carried out in an automatic 
hybridization chamber (Discovery, Ventana Medical sys- 
tem). Images of scanned slides (GenePix 4000A Scan- 
ner- Axon Instruments) were analyzed, spots delimitated 
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and hybridization signals were quantified and trans- 
formed into numerical values by GenePixPro v.3.01 soft- 
ware (Axon). Background noise was rather 
homogeneously distributed and only a few spots were 
saturated at 75%, mainly those corresponding to rRNA. 

Statistical analysis of the data was conducted with the 
R Package Anapuce 2.1 by J. Aubert http://www.agro- 
paristech.fr/ mia/ doku.php?id=productions:logiciels. Nor- 
malization rested on a global lowess regression followed 
by a block correction, after filtering out spots with a sig- 
nal to noise ratio < 3 (including empty spots). Back- 
ground was not subtracted. Differential analysis was 
performed on average values for the triplicate spots 
obtained by the MeanBySpot function. Three models of 
variance were applied: one variance by gene, a common 
variance for all the genes and clusters of genes with 
equal variance (varmixt). Two different multiple testing 
corrections were used to adjust raw P-values, Bonferroni 
correction (which is the most stringent) and False Dis- 
covery Rate of Benjamini and Hochberg, with a nominal 
type I error rate set to 0.05. 

Microarray accession numbers 

The microarray data have been deposited in the Array 
Express database http://www.ebi.ac.uk/arrayexpress/ 
under the accession numbers A-MEXP-2068 (array 
design) and E-MEXP-3238 (experiment). 

Real-time qPCR for quantitation of steady-states 
transcripts 

The mRNAs corresponding to the genes of interest were 
measured by qPCR using SYBR Green fluorescence, 
appropriate specific primers (see additional file 4: list of 
primers) and total first-strand cDNA as template. Con- 
taminating DNA was first eliminated from RNA samples 
using TurboDNA-free from Ambion. Two series of 
reverse transcription and qPGR experiments were car- 
ried out in this study, for monitoring temporal expres- 
sion and overexpression of sigH on the one hand and 
for validation of microarray data on the other hand. 
Total first strand cDNA was produced with random 
hexamer primers (Random Primer 6 5'd(N6)3', Biolabs) 
using either PowerScript Reverse Transcriptase (Glone- 
tech) or PrimeScript Reverse Transcriptase (Takara). 
The quality of each template cDNA was checked using 
the Bioanalyzer 2100 (Agilent). qPGR was performed 
using specific primers (75-100 nM each) according to 
the recommended protocol for each SYBR Green mix 
used (SYBR Green MasterMix 2X from ABgene or 
MESA GREEN MasterMix from Eurogentec). Reactions 
were run on an ABI PRISM 7900 HT instrument 
(Applied Biosystems) or a Mastercycler Realplex 2 S 
instrument (Eppendorf) using 40 cycles of denaturation 
at 95°C for 15 s and extension at 60°C for 1 min. The 



cycles were preceded by DNA polymerase activation at 
95°C and followed by a denaturation cycle to check the 
specificity of the PGR products. Mean Ct obtained for 
studied genes were between 16 and 28.5, with the 
exception of comC and dprA in WT strain at 31 and 
32.9 respectively (in the same time 'No Template Con- 
trols' gave no signal after 34 cycles). 

Primers were designed with Primer Express 2 (Applied 
Biosystems) or Primer 3 http://frodo.wi.mit.edu/primer3 
and validated by determining slopes of standard curves 
for PGR efficiencies between 90% and 100%. In this con- 
text, we used the 2' ' method to express results as 
fold change in the expression of each gene of interest 
relative to a calibrator sample and a reference gene used 
as an internal control for normalization of the results 
[55]. The stability of transcription of the chosen refer- 
ence gene Idh was checked by standard curves per- 
formed for all environmental conditions used in this 
study. Unless otherwise indicated, quantitation experi- 
ments were performed with three independent samples, 
each well being duplicated two or three times. Values 
are expressed as mean ± standard deviation. 

Viability and UV assays 

Viable bacteria were counted by plating serial dilutions 
on MRS agar and incubating at 30°G for one to four 
days. For mixed cultures, classical enumeration on MRS 
supplemented with Xgal (5-bromo-4-chloro-3-indolyl-P- 
D-galactopyranoside, 0.04 g.l '^) distinguished sigH{hy)* 
(white) from sigH(wtY (blue) as well as sigH{nul) (white) 
from 23 K lacLM + (blue). For other tests, sampling for 
stationary phase survival in MGD was done after 6-8 
hour culturing which corresponds to growth arrest, then 
once or twice a day. In these cases, comparative enu- 
meration was performed by depositing drops (5 [il) of 
serial decimal dilutions for each strain on an agar plate. 
UV resistance was examined by exposing bacteria freshly 
plated on MRS medium to 254 nm UV-light (VL-15 G, 
Apelex) with fluences of 40 to 120 J/m^ (by step of 20) 
measured by the radiometer VLX-3 W equipped with a 
254 nm sensor (Vilber Lourmat, France). For that pur- 
pose, drops (5 fil) of serial dilutions of bacterial suspen- 
sions were rapidly plated after growth in MGD medium, 
one hour after addition of 30 |jM GUSO4 or at the same 
time without addition of GUSO4. 

Natural genetic transformation 

Exogenous DNA used in this study comes from plasmid 
(pVI1056, pRV620, and pGKV259) and L. sakei chromo- 
some (strain RV2000), and confers resistance to chlor- 
amphenicol (10 mg.r^) to recipient bacterium. RV2000 
is a derivative of L. sakei 23 K in which the cflfl94 gene 
interrupts Idh [56]. pVI1056 (Van de Guchte in [27]) is 
a 7.5 kb pIPSOl-derived shuttle plasmid (theta- 
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replicating), known to be replicative in L. sakei. pRV620 
is a 5.6 kb shuttle plasmid derived from theta-replicating 
plasmid pRV500 of L. sakei [27]. pGKV259 is a broad 
host range 5 kb rolling circle plasmid [57]. Plasmids 
were purified from E. coli TGI using Qiagen Plasmid 
preparation Kits, and checked by electrophoresis on 
agarose gel for the presence of multimers. 10 ng of plas- 
mids pVI1056 and pRV620 were reportedly able to 
transform B. subtilis naturally competent cells [27]. For 
transformation tests with L. sakei, sigH{hyY overexpres- 
sion strain was cultivated in MCD as described above. 
After 30 to 60 min induction with 30 [iM CUSO4 (at 
usual ODeoo of 0.4 and at ODgoo of 0.2 and 0.9 when 
indicated), aliquots of 100 \A of cell suspension were 
mixed with 100 ng pVI1056 DNA in a microtube and 
incubated for one hour at 30°C. Suspensions were then 
plated on selective MRS medium and incubated for sev- 
eral days at 30°C. As sigHiYiyT strain already contained 
a plasmid, its transformability with incoming plasmid 
was verified by electroporation. Transformation tests on 
plates with other L. sakei strains were done as follows. 
23 K, 64 K [plasmid-cured 64], 332 F [pRV500-cured 
332], 160 K and LTH675 [20,52,58] were cuhivated in 
liquid MRS and MCD medium until late exponential 
phase and plated on the same solid medium supplemen- 
ted with 10 mg.r^ chloramphenicol. Drops of pRV620 
and pGKV259 (50 ng each), and RV2000 chromosome 
(500 ng) were deposited on the plates which were then 
incubated at the indicated temperatures. 

Additional material 



Additional file 1: Alignment of four a"-group sigma factors 

Additional file 2: Genotype of L. sakei strains affected in sigH 

Additional file 3: Competence DNA uptake machinery of 6. subtilis 
and comparison with L. sakei 

Additional file 4: List of primers 



Acknowledgements 

ICE core facilities of INRA at Jouy-en-Josas (J.-C. Helbling, S. Makhzami, E. 
Rebours and P. Martin) and the Biochips platform of Toulouse-Genopole (V. 
Le Berre and collaborators) are acknowledged for their contribution to this 
work. E. Le Moel, J. Elor, and A. Le Neve are acknowledged for participating 
in the study during their training period. We thank our colleagues of Micalis, 
from ELEC team (F. Baraige, S. Chaillou, M.-C. Champomier Verges, M. Daty, 
A. Goubet and I. Lucquin) and other teams (C. Gautier, E. Guedon, I. 
Guillouard and M.-A. Petit) for providing material or advice in experiments; S. 
Chaillou and M. Zagorec for providing the L sakei oligoset; C. Delorme and 
C. Neuveglise for help with polymorphism analysis and phylogeny; M. Van 
de Guchte for help with English. Authors are grateful to J.-P. Claverys for 
helpful discussions on genetic competence, and indebted to A. Gruss for 
improvement of the manuscript. This work was supported by INRA funding. 

Author details 

'UIV1R1319 Micalis, INRA E-78350, Jouy-en-Josas, Erance. ^UMR Micalis, 
AgroParisTech, INRA E-78350, Jouy-en-Josas, Erance. ^UMR1313 Genetique 
Animate et Biologie Integrative, plateforme ICE, INRA E-78350, Jouy-en-Josas, 



Erance. ^Conceptus SAS, 50 avenue de Saint Cloud, F-78000 Versailles, 
Erance. 

Authors' contributions 

SS participated in the design of the study, participated in the sequence 
alignments, carried out construction and characterization of overexpression 
strain and carried out part of the qPCR analysis. CB participated in the 
design of the qPCR analysis. AMCLC conceived and participated in the 
design of the study, carried out and supervised the rest experiments, and 
wrote the manuscript. All authors read and approved the final manuscript. 

Received: 8 June 201 1 Accepted: 12 March 2012 
Published: 12 March 2012 

References 

1. GruberTM, Gross CA: Multiple sigma subunits and the partitioning of 
bacterial transcription space. Annu Rev Microbiol 2003, 57:441-466. 

2. Staron A, Sofia HJ, Dietrich S, Ulrich LE, Liesegang H, Mascher T: The third 
pillar of bacterial signal transduction: classification of the 
extracytoplasmic function (ECF) sigma factor protein family. Mol IVIicrobiol 
2009, 74(3):557-581. 

3. Lonetto M, Gribskov M, Gross CA: The sigma 70 family: sequence 
conservation and evolutionary relationships. J Bacteriol 1 992, 
174(12)3843-3849. 

4. Paget MS, Helmann JD: The sigma70 family of sigma factors. Genome Biol 
2003, 4(1):203. 

5. Britton RA, Eichenberger P, Gonzalez-Pastor JE, Fawcett P, Monson R, 
Losick R, Grossman AD: Genome-wide analysis of the stationary-phase 
sigma factor (sigma-H) regulon of Bacillus subtilis. J Bacteriol 2002, 
184{17):4881-4890 

6. Hilbert DW, Piggot PJ: Compartmentalization of gene expression during 
Bacillus subtilis spore formation. Microbiol IVlol Biol Rev 2004, 68(2):234-262. 

7. Grossman AD: Genetic networks controlling the initiation of sporulation 
and the development of genetic competence in Bacillus subtilis. Annu 
Rev Genet 1995, 29477-508. 

8. Lazazzera BA, Kurt:ser IG, McQuade RS, Grossman AD: An autoregulatory 
circuit affecting peptide signaling in Bacillus subtilis. J Bacteriol 1999, 
181(17):5193-520O 

9. Albano M, Hahn J, Dubnau D: Expression of competence genes in 
Bacillus subtilis. J fiacferio/ 1 987, 169(7):31 10-31 1 7. 

10. Schultz D, Wolynes PG, Ben Jacob F, Onuchic JN: Deciding fate in adverse 
times: sporulation and competence in Bacillus subtilis. Proc Natl Acad Sci 
USA 2009, 106(50):21027-21034. 

1 1 . Nies DH: Incidence and function of sigma factors in Ralstonia 
metallidurans and other bacteria. Arch Microbiol 2004, 1 8(4):255-268. 

12. Morikawa K, Inose Y, Okamura H, IVlaruyama A, Hayashi H, Takeyasu K, 
Ohta T: A new staphylococcal sigma factor in the conserved gene 
cassette: functional significance and implication for the evolutionary 
processes. Genes Cells 2003, 8(8):699-712. 

13. Claverys JP, Mart:in B: Bacterial "competence" genes: signatures of active 
transformation, or only remnants? Trends iVIicrobiol 2003, 1 1(4):161-165. 

14. Kovacs AT, Smits WK, Mironczuk AlVl, Kuipers OP: Ubiquitous late 
competence genes in Bacillus species indicate the presence of 
functional DNA uptake machineries. Environ Microbiol 2009, 
11(8):1 91 1-1922. 

15. Luo P, Morrison DA: Transient association of an alternative sigma factor, 
ComX, with RNA polymerase during the period of competence for 
genetic transformation in Streptococcus pneumoniae. J Bacteriol 2003, 
185(1):349-358. 

16. Chaillou S, Champomler-Verges MC, Cornet M, Crutz-Le Coq AM, 
Dudez AM, IVlartIn V, Beaufils S, Darbon-Rongere E, Bossy R, Loux V, et ah 
The complete genome sequence of the meat-borne lactic acid 
bacterium Lactobacillus sakei 23 K. Nat Biotechnoi 2005, 23(1 2):1 527-1 533. 

17. Jones RJ, Wiklund E, Zagorec M, Tagg JR: Evaluation of stored lamb bio- 
preserved using a three-strain cocktail of Lactobacillus sakei. IVleat Sci 
2010 86(4):955-959. 

18. Vermeiren U Devlieghere E, Debevere J: Evaluation of meat born lactic 
acid bacteria as protective cultures for the biopreservation of cooked 
meat products. Int J Food l<Aicrobioi 2004, 96(2):149-164. 

19. Plewniak E, Jeanmougin E, Higgins DG, Thompson JD, Gibson TJ: The 
CLU5TAL X windows interface: flexible strategies for multiple sequence 



Schmid et al. BMC Microbiology 2012, 12:32 Page 14 of 14 

http://www.biomedcentral.eom/1 471 -2 1 80/1 2/32 



alignment aided by quality analysis tools. Nucleic Acids Res 1 997, 
25(24):4876-4882. 

20. Chaillou S, Daty M, Baraige F, Dudez AM, Anglade P, Jones R, Alpert CA, 
Champomier-Verges MC, Zagorec M: Intraspecies genomic diversity and 
natural population structure of the meat-borne lactic acid bacterium 
Lactobacillus sakei. AppI Environ Microbiol 2009, 75(4):970-980. 

21. Wydau S, Dervyn R, Anba J, Dusko Ehrlich S, Maguin E: Conservation of 
key elements of natural competence in Lactococcus lactis ssp. FEMS 
Microbiol Lett 2006, 257(1 ):32-42. 

22. Morikawa K, Ohniwa RL, Kumano M, Okamura H, Saito S, Ohta T: The sigH 
gene sequence can subspeciate staphylococci. Diagn Microbiol Infect Dis 
2008, 61(4):373-380. 

23. Stentz R, Loizel C, Malleret C, Zagorec M: Development of genetic tools 
for Lactobacillus sakel: disruption of the B-galactosidase gene and use of 
lacZ as a reporter gene to study regulation of the putative copper 
ATPase, AtkB. AppI Environ Microbiol 2000, 66(10):4272-4278. 

24 Weir J, Predich M, Dubnau E, Nair G, Smith I: Regulation of spoOH, a gene 
coding for the Bacillus subtllls sigma H factor. J Bacteriol 1991, 
173(2)521-529. 

25. Desroche N, Beltramo C, Guzzo J: Determination of an internal control to 
apply reverse transcription quantitative PCR to study stress response in 
the lactic acid bacterium Oenococcus oenl. J Microbiol Metliods 2005, 
60(3):325-333. 

26. Tao L, Wu X, Sun B: Alternative sigma factor sigmaH modulates 
prophage integration and excision in Staphylococcus aureus. PLoS Pathog 
2010, 6(5):el 000888. 

27. Crutz-Le Coq AM, Zagorec M: Vectors for lactobacilli and other Gram- 
positive bacteria based on the minimal replicon of pRVSOO from 
Lactobacillus sakei. Plasmid 2008, 60(3):2 12-220. 

28. Chen I, Dubnau D: DNA uptake during bacterial transformation. Nat Rev 
Microbiol 2004 2(3)241-249. 

29. Haldenwang WG: The sigma factors of Bacillus subtilis. Microbiol Rev 1 995, 
59(1):l-3a 

30. Blomqvist T, Steinmoen H, Havarstein IS: Natural genetic transformation: 
A novel tool for efficient genetic engineering of the dairy bacterium 
Streptococcus therrDophilus. AppI Environ Microbiol 2006, 72(1 0):675 1-5756. 

31. Biornstad TJ, Havarstein LS: CIpC acts as a negative regulator of 
competence in Streptococcus thermophilus. Microbiology 201 1, 157(Pt 
6):1 676-1 684. 

32. Berka RM, Hahn J, Albano M, Draskovic I, Persuh M, Cui X, Sloma A, 
Widner W, Dubnau D: Microarray analysis of the Bacillus subtilis K-state: 
genome-wide expression changes dependent on ComK. Mol Microbiol 
2002 43(5):1331-1345. 

33. Mortier-Barriere I, Velten M, Dupaigne P, Mirouze N, Pietrement 0, 
McGovern S, Fichant G, Martin B, Noirot P, Le Cam E, et al: A key 
presynaptic role in transformation for a w/idespread bacterial protein: 
DprA conveys incoming ssDNA to RecA. Cell 2007, 130(5):824-836. 

34 Brown NL, Stoyanov JV, Kidd SP, Hobman JL: The MerR family of 
transcriptional regulators. FEMS Microbiol Rev 2003, 27(2-3):145-163. 

35. Peterson S, Cline RT, Tettelin H, Sharov V, Morrison DA: Gene expression 
analysis of the Streptococcus pr^eumor^iae competence regulons by use 
of DNA microarrays. J Bacteriol 2000, 182(21):6192-6202. 

36. Claverys JP, Martin B, Polard P: The genetic transformation machinery: 
composition, localization, and mechanism. FEMS Microbiol Rev 2009, 
33(3):643-656. 

37. Lindner C, Nijiand R, van Hartskamp M, Bron S, Hamoen LW, Kuipers OP: 
Differential expression of two paralogous genes of Bacillus subtllls 
encoding single-stranded DNA binding protein. J Bacteriol 2004, 
186(4):1 097-1 105. 

38. Cardan R, Besset C, Guillot A, Gitton C, Monnet V: The oligopeptide 
transport system is essential for the development of natural 
competence in Streptococcus therrr]ophilus strain LMD-9. J Bacteriol 2009, 
191(14):4647-4655. 

39. Maughan H, Redfield RJ: Tracing the evolution of competence in 
Haemophilus Influenzae. PLoS ONE 2009, 4(6):e5854. 

40. Johnsborg 0, Eldholm V, Havarstein LS: Natural genetic transformation: 
prevalence, mechanisms and function. Res Microbiol 2007, 
158(10):767-778. 

41. Hofer F: Transfer of lactose-fermenting ability in Lactobacillus lactis. New 

Zealand Journal of Dairy Science and Technology 1985, 20(3):179-183. 



42. Meibom KL, Blokesch M, Dolganov NA Wu CY, Schoolnik GK: Chitin 
induces natural competence in Vibrio cholerae. Science 2005, 
310(5755):1824-1827. 

43. Shaw AJ, Hogsett DA, Lynd LR: Natural competence in 
Thermoanaerobacter and Thermoanaerobacterium species. AppI Environ 
Microbiol 2010, 76(14):4713-4719. 

44. Nijiand R, Burgess JG, Errington J, Veening JW: Transformation of 
environmental Bacillus subtilis isolates by transiently inducing genetic 
competence. PLoS ONE 2010, 5(3):e9724. 

45. Woodbury RL, Wang X, Moran CP Jn Sigma X induces competence gene 
expression in Streptococcus pyogenes. Res Microbiol 2006, 157(9):851-856. 

46. Mashburn-Warren L, Morrison DA, Federle MJ: A novel double-tryptophan 
peptide pheromone controls competence in Streptococcus spp. via an 
Rgg regulator. Mol Microbiol 2010, 78{3):589-606. 

47. Metzger Z, Dotan M, Better H, Abramovitz I: Sensitivity of oral bacteria to 
254 nm ultraviolet light. Int Endod J 2007, 40(2):120-1 27. 

48. Phillips ZE, Strauch MA: Bacillus subtilis sporulation and stationary phase 
gene expression. Cell Mol Life Sci 2002, 59(3):392-402. 

49. De Man JC, Rogosa M, Sharpe ME: A medium for the cultivation of 
lactobacilli. J AppI Bacteriol 1960, 23:130-135. 

50. Lauret R, Morel-Deville F, Berthier F, Champomier-Verges M, Postma P, 
Ehrlich SD, Zagorec M: Carbohydrate utilization in Lactobacillus sake. AppI 
Environ Microbiol 1996, 62{6):1 922-1 927. 

51. Hungate RE: A roll tube method for the cultivation of strict anaerobes. In 
Methods in Microbiology Volume 3B. Edited by: Norris JR, Robbons DW. 
London: Academic Press; 1969:117-132. 

52. Alpert CA, Crutz-Le Coq AM, Malleret C, Zagorec M: Characterization of a 
theta-type plasmid from Lactobacillus sakei: a potential basis for low- 
copy-number vectors in lactobacilli. AppI Environ Microbiol 2003, 
69{9):5574-5584 

53. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary 
Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 
24{8):1 596-1 599. 

54. Herve-Jimenez L, Cuillouard I, Guedon E, Gautier C, Boudebbouze S, Hols P, 
Monnet V, Rul F, Maguin E: Physiology of Streptococcus thermophilus 
during the late stage of milk fermentation with special regard to sulfur 
amino-acid metabolism. Proteomics 2008, 8(20):4273-4286. 

55. Livak KJ, Schmittgen TD: Analysis of Relative Gene Expression Data Using 
Real-Time Quantitative PCR and the 2-AACT Method. Methods 2001, 
25{4):402-408. 

56. Malleret C, Lauret R, Ehrlich SD, Morel-Deville F, Zagorec M: Disruption of 
the sole Idh gene in Lactobacillus sakei prevents the production of both 
L- and D-lactate. Microbiology 1998, 144:3327-3333. 

57. Kok J: Special-purpose vectors for lactococci. In Genetics and molecular 
biology of streptococci, lactococci and enterococci Edited by: Dunny GM, 
Cleary PP, McKay LL. Washington, D.C: American Society for Microbiology; 
1991:97-102. 

58. Berthier F, Zagorec M, Champomier-Verges M, Ehrlich SD, Morel-Deville F: 
High-frequency transformation of Lactobacillus sake by electroporation. 

Microbiology 1996, 142:1273-1279. 



doi:10.1 186/1471-2180-12-32 

Cite this article as: Schmid et al: Alternative sigma factor a" activates 
competence gene expression in Lactobacillus sakei. BMC Microbiology 
2012 12:32 



Submit your next manuscript to BiolVled Central 
and tal<e full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



BioMed Central 



